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The gauge problem in the so-called strong-field approximation (SFA) describing atomic or molec- 
ular systems exposed to intense laser fields is investigated. Introducing a generalized gauge and 
partitioning of the Hamiltonian it is demonstrated that the S-matrix expansion obtained in the 
SFA depends on both gauge and partitioning in such a way that two gauges always yield the same 
S'-matrrx expansion, if the partitioning is properly chosen. 

PACS numbers: 32.80.Rm, 33.80.Rv 

I. INTRODUCTION 

Gauge-invariance is one of the fundamental concepts of electrodynamics. As a consequence it is, e.g., possible to 
formulate the interaction of charged particles with electromagnetic fields in different gauges. Although the choice of 
the gauge clearly influences parameters like the scalar or the vector potential, all physical quantities (observables) 
are independent of the gauge, if a complete treatment is performed. On the other hand, an approximate treatment 
often leads to gauge-dependent predictions for physical observables. One prominent example is given by the so-called 
strong-field approximation (SFA) for describing atomic or molecular systems exposed to intense laser fields that is 
also known as Keldysh-Faisal-Reiss (KFR) theory [U, 0, Q . It is based on an (infinite) series expansion of the exact 
S matrix describing the interaction of an atomic system with a laser pulse. Besides a phase factor, the S matrix 
obtained from a calculation of the complete expansion should thus be independent on the chosen gauge, provided 
the series converge. Consequently, the corresponding observable transition probabilities (obtained from the squared 
absolute values of the S'-matrix elements) should be gauge independent. 

This gauge-independence of physical observables is usually lost, if only a truncated series is considered. This is the 
case for the SFA which is defined as the first term of the ^-matrix expansion. In a number of investigations it has 
been shown that transition probabilities or rates predicted by either the length- or the velocity-gauge formulation of 
the SFA differ easily by one or two orders of magnitude for experimentally relevant laser parameters . Furthermore, 
it was shown that the velocity-gauge SFA rate does not converge to the tunneling limit for weak fields, if long- 
range Coulomb interactions are present @ . Recently, it was also demonstrated that there are pronounced qualitative 
differences between the energy distributions of the electrons ejected from, e. g., the 2p state of hydrogen atoms exposed 
to intense laser fields, if they are calculated within the SFA and either the length or the velocity gauge |6|]. The recent 
extensions of the SFA to molecular systems in velocity gauge , length gauge [1] , or dressed length gauge @ indicate 
similar or even more pronounced gauge dependencies for molecules. These findings have intensified the discussions 
whether the formulation of the SFA in one of the two gauges may be "superior" to the other. One approach to answer 
this question is rather pragmatic and based on a direct comparison of the SFA predictions in both gauges to either 
exact (numerical) solutions of the full time-dependent Schrodinger equation [13] or to experimental results [ill, ■ 
Clearly, if there is no a priori reason that one of the two SFA formulations is superior to the other, the conclusions 
may vary depending on the considered atomic or molecular system (even its quantum state) as well as on laser-pulse 
parameters. 

A second line of argumentation in favor of one of the two gauges is based on the question of "universality" . 
For example, it has been argued that an evident limitation of the length-gauge formulation is the fact that in this 
case the predicted observables depend only on the scalar potential and thus a one-dimensional parameter, while 
the full description of an electromagnetic field requires in principle more than one dimension [l3|. Very recently, 
Faisal proposes a "gauge-invariant" intense-field iS-matrix theory that yields equal transition probabilities in length 
or velocity gauge; independent of the order of expansion [l3|. Consequently, Faisal claims that his theory overcomes 
the above-mentioned long-standing discrepancy between the SFA in the two gauges. According to the findings in [l3| 
the "gauge-invariant" and thus universal S'-matrix theory appears to be equivalent to the traditional length-gauge 
formulation. Specifically, the first-order term reproduces exactly the Keldysh result '1] which was obtained in length 
gauge. In view of the popularity of the SFA for describing atomic and molecular ionization in intense laser fields (see, 
e.g., (l5| or and references therein), this is an important result. 

Besides the evident appeal of a universal S'-matrix formulation which would provide an end to the long-lasting 
debates on the choice of the appropriate gauge (for a very recent example, see [l3, [3), it appears, however, quite 
surprising that such a formulation should exist. In fact, Faisal derives in [19] an alternative S-matrix expansion that 
reproduces in length and velocity gauge the traditional velocity-gauge result. This alternative series expansion is 
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obtained by a different partioning of the Hamiltonian. The particular results obtained in fl^ and fl^ immediately 
lead to the question of their generality. In the present work it is shown that it is possible to achieve an S'-matrix 
expansion in agreement to any traditional SFA formulation in either length, velocity, or radiation gauge for an 
arbitrary choice of the gauge, if the Hamiltonian is correspondingly partitioned. In fact, introducing a generalized 
gauge transformation that includes the mentioned particular gauges as special cases, an in principle infinite set of 
different S'-matrix expansions can be formulated. All of them can be shown to be achievable within any gauge as long 
as a proper partioning of the Hamiltonian is performed. Furthermore, the introduction of the generalized gauge allows 
to clearly demonstrate, how the choice of the gauge and the partitioning of the Hamiltonian describing the atomic or 
molecular system exposed to a laser field are connected with each other. This provides a much deeper insight in the 
gauge-problem of SFA that in fact turns out to be more properly described as an expansion problem. 

In order to provide a clear definition of terms and notations, the following section gives a brief overview over the 
gauges relevant to the present work, including brief discussions on (local) gauge invariance, form invariance, and 
the dipole approximation. Most importantly, a generalized gauge is introduced in Sec. Ill Gl Sec. [Hi] discusses the 
wavefunction of a free-electron in an electromagnetic field and the influence of the chosen gauge on it. Equipped 
with these prerequisites, Sec. lIVI discusses the S'-matrix theory in different gauges. It represents thus the main results 
of this work in which it is demonstrated how various results (S'-matrix expansions) can be obtained using different 
combinations of gauge and partitioning. 



II. GAUGES 



In the validity regime of the SFA the laser intensities are so high that the photon density is also very high. Under 
these circumstances the number of photons can be treated as a continuous variable and the field can be described 
classically by using Maxwell's equations. Therefore, in the following a semi-classical theory will be used in which the 
radiation field is treated classically, but the atomic or molecular system is described using quantum mechanics. The 
influence of the quantum system on the external field is also neglected. 



A. Local gauge invariance 

A classical electromagnetic field is described by electric and magnetic field vectors, F{r,t) and B(r, i) or, alter- 
natively, by the scalar and vector potentials $(r,t) and A{r,t). Consider a system consisting of an electron in a 
electrostatic potential U{r) created by a nucleus (or some nuclei) which interacts with an external electromagnetic 
field. In semi-classical theory, the evolution of the system is governed by the time-dependent Schrodinger equation that 
satisfies local gauge invariance and is given in the coordinate representation with the minimal-coupling Hamiltonian 

by 



^^*x(r-0 = ^ o Pc + Ax(r,i) ~ <^^{r,t) + U{r) } ^^{r,t) (1) 



dt ' ' [2 

where the subscript denotes the used gauge x ^^'^ the operator of canonical momentum is given independently of 
the gauge as p^. = — iV; a consequence of the definition of the minimal-coupling Hamiltonian. The word invariance 
means that if the wavefunction and both potentials and are simultaneously transformed into a new gauge 
x' using the transformation recipes 

= ^xlr^Oet'^^^^'*""'*^] *x(i",i) = '^^'{r,t)e\-''^x'^>c^''^'^^ (2) 

A^,(r,i)=A^(r,t)-VT^^x'(r,t) A^{r,t) ^ A^,{r,t) - yT^,^^{r,t) (3) 

d d 



equation ([T|) is transformed into itself, but with "^x'' -^x ~^ -^x'' *^x ^ *^x'- '^^'^ transformation functions 

T^^^'(r,i) satisfy the relations 



Of course, physical quantities as the probability P{r,t) or the electric and magnetic fields (F(r, i) and B(r, i)) are 
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jauge independent, i. e. 



P(r,<) = |vl/^(r,t)p = |*^,(r,i)p (6) 



B{r,t)^VxA^{r,t)^VxAy{r,t). (8) 



F(r,t) = -V<i>x(r,i) - ^A^(r,t) = -V$^-(r,i) - ^^^^'(r,*) (7) 



B. Form-invariant physical quantities 



A form- invariant quantity is defined as a quantity whose corresponding operator — G{A^, is form invariant 
under a unitary transformation Tx^x' ~ ^^P [*'^x^x'('': '^)]^ i- S-, 



GiAx',^, 



(9) 



The important difference between physical and non-physical quantities lies in the gauge invariance of the eigenvalues. 
The eigenvalues of a physical quantity (an observable) are identical in all gauges, whereas the eigenvalues of non- 
physical quantities may depend on the chosen gauge. It can be shown that the eigenvalues of a form-invariant quantity 
are gauge- invariant, i. e. a physical quantity must be form-invariant. Consider the operator with the eigenvalues 
gn and the corresponding eigenstates \x,n), 



Gx\x,n) = gn\x,n) 



(10) 



The operator G^' obtained by the unitary transformation ([9]) has the eigenstates = T^^t^'|x, n) and the same 

eigenvalues g„ 



G'x'lx','^) = i^x^x'^x'^ x^x')i^x^x'\x,n)) ^ Tx^x' 9n\x,n) ^ g,,\x' ,n) 



(11) 



Both and <i>^ are not operators that correspond to physical quantities. The canonical momentum is also not 
form- invariant , 
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and, therefore, does not represent a physical measurable quantity, whereas the mechanical momentum 

TTx = Pc + A;<:(r,t) 

is form-invariant. 



Tx^x'^^xT^^x' 



VTx^x' (r, t)+Ax = V, + Ax' = tt. 



(12) 
(13) 

(14) 



and represents an observable. The electrostatic potential U (r) also represents a physical quantity. The instantaneous 
energy operator of the system, f^, 



^x = ^ 
^ 2 



1 r 



Pc + Ax(r,i) +C/(r) 



C/(r) 



(15) 



represents a physical quantity, since it is a function of only physical quantities: tt^ and U{r). 
It is essential that both the total Hamiltonian of the system. 
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Pc + Ax(r, t)J - ^xir, t) + U{v) =£x- t), 
and the field-free Hamiltonian, 

Ho^^ + [/(r), 

are not form- invariant, although they determine the evolution of a physical system. 



(16) 



(17) 
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C. Radiation gauge 



The radiation gauge (labeled in the following by the subscript R) is convenient, if no sources are presented. It is 
defined by the relations 

V ■ Ar = 0, $R = . (18) 

Therefore, one has V • (Ar'S) = Ar(V^') + (V • Ar)^!' — Ar(V\1/) and the total Hamiltonian of the system possesses 
the form 

Hr = p2/2 + Ar ■ p, + A^/2 + [/ . (19) 

The vector potential in this gauge satisfies the wave equation and can be given by a superposition of monochromatic 
plane wave solutions. In the following, the vector function A(r,t) specifies the vector potential AR(r, i) in radiation 
gauge, so that 

F(r,t) = -|-A(r,t), B(r, t) = V x A(r, i). (20) 



D. Dipole approximation 

If the wavelength of the considered radiation is sufficiently long, the spatial variation of the radiation field across 
the system can be neglected. Assuming A(r,t) w A{t) one obtains 

F = Fit) = -^^^, B = 0. (21) 
dt 

There exist two gauges, length and velocity gauge, which are extensively used in the context of the dipole approxima- 
tion. It is, however, known that the conditions for a breakdown of the dipole approximation differ in different gauges. 
In order to explicitly highlight the terms which are neglected in the subsequent use of the dipole approximation, the 
general definitions of length and velocity gauges are given in the following subsections. 



E. Velocity gauge 

The velocity gauge (labeled in the following by the subscript V) is often used to remove the square of the vector 
potential from the total Hamiltonian. Note, the name velocity gauge is frequently used, if in fact the radiation gauge 
is meant. The reason is the equivalence of these two gauges in the weak- field limit in which the term can be 
neglected (with respect to the remaining terms). For strong fields this term is, however, not negligible and it is thus 
important to distinguish between these two gauges. The velocity gauge can in fact be obtained from the radiation 
gauge by means of the transformation function 

Tn^yir,t)^^j\\r,t')dt' = P{r,t) . (22) 

The vector and scalar potentials are given in the velocity gauge by 

Av = A - V/3, $v = AV2 (23) 
and the total Hamiltonian can be rewritten as 

Hv = Pe/2-l-AvPc + iA/3 + [/ (24) 

where the identity V • Ay = —A/3 has been used. 

Within the dipole approximation, the vector potential has the simpler form Ay = A(i) and Eq. ((24|) reduces to 



Hv - pI/2 + U + Ait) ■ p. 



(25) 



5 



F. Length gauge 

Another useful and popular gauge is the length gauge (labeled in the following by the subscript L). It is obtained 
from the radiation gauge, if the transformation function 

TR,^L(r,t) = A(r,t).r (26) 

is applied. This leads to the length-gauge vector and scalar potentials 

Al = A- V(A-r) = -iL X A, $l = -F • r (27) 

where L = r x p^. is the canonical angular momentum. Within the dipole approximation, the vector potential in 
length gauge is simply zero. This leads to a simple form of the total Hamiltonian that is given in length gauge by 

HL = PeV2 + C/ + F(i)-r. (28) 

Note, that in this case the operators of the mechanical and canonical momentum coincide, ttl = Pc. 

G. Generalized gauge 

It can be shown that both the length and velocity gauges are only particular cases of a generalized gauge defined 
by an arbitrary set of parameters, X — {a;i,X2}. This gauge, which will be referred to as X gauge, is obtained via 
the transformation 

TR^x(r,i) = a;iA(i) •r + X2/3(t) (29) 

In the X gauge the vector and scalar potentials are given as 

Ax(t) - (l-a;i)A(t), $x = -a;iF(i)T + a;2A2(t)/2 . (30) 

This leads to the total Hamiltonian 

iix = Pc/2 + U+{1- xi)A(t) • p, + XI Fit) • r + [(1 - xif - X2] A^{t)/2. (31) 

In particular, X = {0, 0}, X — {0, 1}, and X — {1, 0} correspond to the total Hamiltonians in radiation gauge (|19|1. 
velocity gauge ([25| . and length gauge ((28|) . respectively. 

III. FREE ELECTRON IN A LASER FIELD 
A. Volkov wave function 

Consider now a free electron in the presence of a laser field described by A{t). In velocity gauge the TDSE reads 

j^*v(r, t) = [pI/2 + Ait) ■ pj ^^(r, t). (32) 

The solution of Eq. (j32p is the (non-relativistic) Volkov wavefunction 

5'v,k(r, t) = exp ik • {r - a{t)} - iEkt , Ek ^ k'^/2. (33) 

where a{t) — J* A{t') dt' and (5-function normalization is used. It is important to emphasize the physical meaning of 
the vector k and the scalar Ek that are usually referred to as electron momentum and electron energy, respectively. 
In fact, k is the mean value of canonical momentum in velocity gauge, (Pc)v = k. The mean value of mechanical 
momentum is on the other hand time-dependent and equal to k -I- A(t). As was discussed above, this latter value is 
gauge-independent and will be denoted in the following by tt; thus one has tt = 7r(k, <) = k + A(t). To stress the 
independence on a gauge, the vector function tt has no subscript, in contrast to the operators ttr, ttv, or ttl- The 
mean value f (k, t) — 7r^(k, t)/2 of the instantaneous energy operator (fT5|) is thus also time-dependent. 
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Using the identity 



1 /■* 

- / £:(k, <') dt' ^Ekt + k- a{t) + fi{t) 



the Volkov wavefunction can be written as either 



(34) 



^R,k(r, = exp ik • r - i / £{\l, t') dt 



(35) 



*L.k(r, t) = exp i7r(k, t)-v-i / £(k, t') dt 



in radiation and length gauge, respectively. 
In X gauge one finds 



4'jf,k(r, t) = exp ixi A(t) ■ r + ix2/?(t) «'R,k(r, t) = e-'©^(*)e*-- 



(36) 



(37) 



with 



Qxit) = Ekt + k- a{t) - xiA{t) ■r-{x2- l)P{t) 



(38) 



Note, in length gauge k can be seen purely as a parameter, since also the mean value of the canonical momentum 
(P(,)l is equal to tt. Nevertheless, the physical meaning of k becomes evident in the case of a linear-polarized 
monochromatic electromagnetic field. 



B. Linear-polarized monochromatic field 



For a vector potential given by A{t) — —Aq sin tut the corresponding electric field is F{t) — Fq cos uit with Fq = luAq. 
In this case the mean value of the mechanical momentum is equal to 7r(k, i) = k — Fo/ojsmojt. Therefore, k is the 
value of the cycle-averaged tt. An electron with k = quivers in the field around a single point in space. 

The mean value f (k, t) of the instantaneous energy operator (jisp is time-dependent and given by 

1 / F \ F 

£(k,t) = - k -smut] = Ek -k- —sinuit- Up cos2Ljt + Up. (39) 

2 \ uj J uj 

Therefore, its cycle-averaged value is equal to Ek + Up. The pondcromotive energy Up = Fq/{4uj'^) is thus the 
cycle-averaged instantaneous energy of a quivering electron (k = 0) . 



IV. FORMAL S-MATRIX FORMULATION OF SFA 



The following formulation of the S'-matrix theory describing atomic and molecular systems in intense laser fields 
considers the case of a one-electron system for the sake of simplicity. The generalization to an arbitrary number of 
electrons is, however, straightforward. As a starting point the TDSE formulated in X gauge is considered, 

(^i|-Hx(t))|*^(<))=0. (40) 

The electromagnetic field is absent before and after the pulse, i. e. 

A{t) = F{t) = 0, Hjf(t)=H" for t < and t > t/ (41) 
Here, H'' is the ficld-frce Hamiltonian with eigenvalues Ea and eigenvectors |'0a), 

11° = pI/2 + U, H°|^„) - EM- (42) 
(The index a denotes discrete as well as continuum states and is thus itself either discrete or continuous.) 
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To describe the action of the pulse on the system, complete and orthonormal initial- and final-state basis sets are 
introduced. The initial-state basis set is given by |V'a(ii)) = e~'^°*'|'0a) where the phase factor is introduced for 
convenience. The final-state basis set is given by plane waves with momentum k, again for convenience multiplied by 
a phase factor, and depends both on the adopted gauge and on k, 

|*x,k(i/)) =e-'®-(*^)|k). (43) 

The phase (see Eqs. ^ and (j4T|) ) 

&x{tf) - E,tf + k • a{tf) - {X2 - l)P{tf) (44) 

is r independent but depends on the used gauge, k, and the pulse. Note, that the for reasons of convenience introduced 
phase factors add only constant phases in the transition amplitudes and do not alter transition probabilities. 
The probability amplitude of a transition from an initial state \ipa{ti)) to a final state \^ x,v.{tf)) is given by 

^ka - i(*x,k(i/)|Gx(i/,t»)|^a(i^)) (45) 
where the propagator Gx{t^t') is associated with ^i-x{t) by the inhomogeneous equation 

d 



ilx{t)^Gx{t,t')=S{t-t'). (46) 



To obtain a systematic expansion of the transition amplitudes of interest it is convenient to express the total 
propagator Gx of the system in terms of a partial propagator, defined by a partitioning of the total Hamiltonian. 
The choice of the partitioning is made in such a way that the partial propagator can be expressed analytically, i. e. 
the Schrodinger equation with the corresponding partial Hamiltonian is solvable. 

A first class of Hamiltonians that leads to analytical solutions is the one describing a free electron in the field. As 
was discussed in Sec. IIIH such Hamiltonians are gauge-dependent and their solutions are given by Volkov states. The 
partitioning of Hx using the free-electron Hamiltonian in X gauge, H^, is given by 

Rx^&x + U (47) 
The corresponding propagator can be written analytically using the solutions |5'x.k(0)i 

G'xit,t') = -^e{t-t')J2\'^x,wit)){^xMt% (48) 

k 

where 9(x) is the step function. From Eq. (|48|l follows 

-i{^x,uitf)\G'xitf,t) = {^XMit)\ ioyt<tf. (49) 

Another Hamiltonian that can be used for the partitioning is the field-free Hamiltonian H*'. It is, however, only a 
special case of the class of Hamiltonians, which will be referred to as generalized field-free Hamiltonians and will be 
considered in the following subsection. 



A. Generalized field-free Hamiltonian 



The generalized field-free Hamiltonian 

HO ^ e-(-*)H"e--(-*' - = H° + ^A^ - (V7) • + ^ {Vjf - ^ (50) 

is defined with the aid of an arbitrary function 7(r,t) in such a way that it reduces for 7 = to the field-free 
Hamiltonian H" . Solutions of the TDSE with H° can be expanded in terms of the solutions of the TDSE with H" as 

=e'^('-'*)e-'^°*|Va). (51) 

In general, the function 7(r, t) can be chosen independent of the gauge that is used to formulate the TDSE. Consider 
a particular choice of 7(r,t) parameterized by a set of parameters, A — {Ai, A2}, 

7A(r,i) = AiA(t)-r-|-A2/3(t) . (52) 
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The corresponding generalized field-free Hamiltonian (the subscript A is adopted instead of 7a for the sake of notational 
simplicity) is then given by 

fio = H" - Ai A{t) • p, + Ai F{t) -r+iXl- \2)A^{t)/2. (53) 

Note, for all choices of A the Hamiltonian H° gives an equivalent description of the evolution before and after the 
pulse, since for those times both A{t) and F(i) are equal to zero. Different choices of A yield, however, different 
partial propagators during the pulse that can be written analytically as 

(<, t') = ~ie{t - t') J2 l*A,a(i))(*A,a(i')|- (54) 



In order to express the total propagator Gx in terms of G|^, the total Hamiltonian is partitioned in two parts, 

l" + V" 

where the interaction operator is given by 



H^ = H", + (55) 



Vlxit) = {l-x, + Ai)A(t) • p, + (xi - Ai) F(i) • r 

+ [(1-Xl)2-.X2-A?+A2] A2(i)/2 

It is worth reminding that both sets of parameters, X = {xi,X2} and A = {Ai, A2}, are independent of each other. 
Therefore, the same interaction operator can be obtained for different X gauges, if the A parameters are appropriately 
chosen. It can be shown, for example, that 

^R,{-i,o} = ^v,{-i,i} = ^Mo,o} = m ■ r (57) 
^R,{0,-1} = Vl^o,o} - ^L,{i,i} - Mt) ■ Pc (58) 
^R°{o,o} =^v,{o,i} = ^LV2}=A(t)-p, + A2(t)/2 (59) 

Since |4'A,a(t»)) = li^c^ik)), Eq. ^ yields 

iGUt,U)\i,^{U)) = |*A,a(t)), for t > U. (60) 

B. Matrix elements 

It will now be shown that most of the matrix elements of interest depend at most on the two parameters v = {vi, V2} 
with = 1 + Ai — Xi and 'i;2 = A2 — 0:2. 

Indeed, one finds for different matrix elements the relations 

i^xMimxAt)) = (k|e'^"-W|^„), (61) 
(*x,k(t)|C/|«'A,a(t)) = (k|e'^-°Wt/|^„), (62) 

{^X,a'{t)\Vl^{mxAt)) = e'(^°'^^°^*(^a'|V;!'(i)l^a), (63) 



and 



where 



and 



i^xMitWl.imxAt)) = (k|e'"-°(*)v;°(0|Va) (64) 
V°{t) = Vl.it) + Ai(l - XI + \i)A\t) 



= viA{t) • p, + (1 - vi)F{t) ■r+ivf+ V2)A\t)/2 

f^.ka it) = Qx (t) +jx{r,t)-Eo,t 

= {Ek - E^)t + k • ait) + («i - l)A(i) • r + («2 + l)/3(t) 



(65) 



(66) 
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Finally, the matrix element 

(*x,k'W|C^I*x,kW) =c'(^'='-^^)*+'('''-'')-°(*)(k'|[/|k) (67) 

is independent both on gauge and partitioning. 

As a consequence of these properties of the matrix elements the transition amplitude depends only on w, as is shown 
below. 

C. S-matrix series 

The operator Gx{t,t') can be expanded either in terms of the operator G^{t,t'), 

Gx{t,t')^Gl{t,t') + JdhGx{t,h)V°^^{h)Gl{h,t'), (68) 
or in terms of the operator G^^{t,t'), 

Gx{t,t')^G'x{t,t')+ fdhG'x{t,h)UGx{ti,t'). (69) 



Substitution of Eq. dM]) in ([55]) yields 

Gx{t,t') = G'i{t,t') + ldhG'^{t,h)Vl^{h)G'i{h,t') 

(dt2 AhG'x{tM)UGx{t2M)V^^{ti)Gl(ti,t'). 



(70) 



A further substitution of either ([B5)) or in Eq. fTO]) results in a series expansion of Gxit,t'). Inserting this 
expansion in Eq. (|45|) generates the S'-matrix series for the transition amplitude between the initial state and the final 
state to any desired order, 

oo 

^k. = E4o' (71) 

n=0 

with 

4a =i(*x,k(t/)|G°(t/,iOl^a(i.)), (72) 
4a = i y'dti(vI/^,k(t/)|GS,(t/,ti)P^i;,(ii)G^(ti,tO|Va(iO)- (73) 

Depending on whether ([55)1 or ([55]) is substituted in (|7D|) one obtains either 

4a =i / /'diid^2(vI/x,k(i/)|Gf,(i/,t2)C/G",(i2,ti)y^,,(ti)G"^(ii,i.)|^„(t,)) (74) 

or 

4a =i / /diidt2(*x,k(i/)|Gf,(t/,t2)C/Gf,(t2,ii)V^lA(ii)G!^(ii,^*)lV'a(i^)) (75) 



and so on, where the integration is performed in the range U to tf. 
From Eqs. and (HI]) follows 



4 a = (*X,k(t/)|*A,a(i/)) = e'''-°^*^^Va(k) (76) 



where V'a(k) = (k|V'a) is the Fourier transform of ipa- 
From Eqs. (pOl) and (|15)) follows on the other hand 



4a = i /*'di(vI/x,k(Ol^lA(OI*A,a(0) (77) 
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which — using the identity (|64|) — can be reduced to 

4a = i^* dt(k|e'""-Wv;°Wl^a)- (78) 
In an analogous way, Eq. (fM)) can be transformed using (IM)) as 

= j 'dt2^{^x,v.{t2)\U\^x.o.\t2)) j dti(*A,a'(ii)|V'xa(ii)l*A,a(ti)) (79) 
= ^*'dt2^(k|e'""-'(*^)t/|^„0^*dtie'(^<>'-^")*MV^a'|V;°(ii)l^a) 



or Eq. (|75|) can be transformed using ([481) as 

4a / 'dt2^ V(*X,k(i2)|t/Gff(i2,ti)FiA(tl)|^A,a(il)) 

= ^*'d^2^(vI/^,k(^2)|c/|*x,k'(^2))^V(*A^k'(^l)|v^^,A(^l)l*A,a (80) 



= /" ' <lt2'Y,'i'^'^''~^'''^'^^'^^~^'^'°'^'^H^\U\^) f dti(k'|e'"">'' = (*i)v;°(<i)|V'a). 

Continuing in an analogous manner, it can be shown that S^^J for any order n depends only on v. Therefore, 5ka 
itself depends only on v. In the next subsection some particular cases will be considered explicitly. 

D. Particular cases 

As a first example, consider the case studied in [14]. It is obtained using v — {0, 0} where one has 

K'{t)=F{t)-r (81) 

ri.ka(i) = (Ek ~ E^)t + k • a{t) - A{t) ■ r + /3(t). (82) 
This formulation is achieved using the following partitionings for different gauges: 

A = { — 1, 0} in R-gauge, A = { — 1, 1} in V-gauge, A = {0, 0} in L-gauge. (83) 

Since in L gauge the relation 11° — holds, the gauge-invariant formulation with v ~ {0, 0} reproduces the traditional 
SFA in the L gauge. 

However, in an analogous way, the traditional V-gauge SFA is obtained with v = {1,-1} (cf. [l^), where 

V°{t)^A{t)-p, (84) 
n,i,a{t)^iEk-Eo,)t + k-a{t) . (85) 

It can be achieved using the following partitionings for different gauges: 

A — {0, —1} in R-gauge, A ~ {0, 0} in V-gauge, A = {1, —1} in L-gauge. (86) 

In a similar way, the traditional R-gauge SFA is obtained with v = {1,0}, where 

V°{t)=A{t)-p,+A\t)/2 (87) 
ri.ka(i) = {Ek - Ea,)t + k • a{t) + p{t). (88) 

It can be achieved using the following partitionings for different gauges: 

A = {0, 0} in R-gauge, A = {0, 1} in V-gauge, A = {1, 0} in L-gauge. (89) 

Clearly, every S'-matrix expansion (SFA formulation) in one of the "conventional" (length, velocity, or radiation) 
gauges can be obtained by adopting any of the possible gauges, if the partitioning of the Hamiltonian is chosen 
accordingly. 
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V. CONCLUSION 

In this work it is shown how the (infinite-order) iS-inatrix expansion describing atomic or molecular systems in 
intense laser fields depends on the choice of both the gauge and the partitioning of the Hamiltonian. For this purpose 
a generalized gauge as well as a generalized partioning scheme is introduced. They are defined by 4 independent 
parameters {xi and X2 for the gauge as well as Ai and A2 for the partitioning). However, the 5- matrix expansion is 
then shown to depend only on 2 parameters, Wi = 1 + Ai — Xi and V2 = X2 — X2- Clearly, every possible combination 
of the parameters defining the gauge and the partitioning that conserves the values Ai — Xi and A2 — X2 leads to an 
identical 5'- matrix expansion (up to an overall phase factor that cancels when calculating physical observables). 

The present analysis thus shows that one has to be very careful not to consider solely the gauge, since a suitable 
choice of the partitioning may lead to an S'-matrix expansion that is identical to the one obtained in some other 
gauge, as was also demonstrated in [l^ where the velocity-gauge expansion appeared to agree with what is usually 
known as length-gauge S'-matrix expansion. In fact, it is much more appropriate to discuss different expansions than 
different gauges. This avoids for example the rather contradictory terminology of a "gauge-invariant (first-order) KFR 
approximation in the velocity gauge" used in [T3 |. 

Furthermore, the present result demonstrates that there remains an in principle infinite set of S'-matrix expansions 
(characterized by different values of vi,V2) that are only shown to provide the same transition probabilities in the 
limit of an infinite series expansion, if the latter converges. Truncated series like, e.g., the 0th, 1st, or 2nd order 
expansions will, however, in general disagree. Thus the question of a "proper" choice of the expansion in the case 
of truncation remains and can only be clarified by a comparison to either experiment or gauge-independent theory 
(like full solutions of the time-dependent Schrodinger equation). It should be emphasized, however, that there is no 
a priori reason to believe that one expansion is necessarily advantageous to the others. 
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